Entanglement of a two-particle Gaussian state interacting with a heat bath 
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The effect of a thermal reservoir is investigated on a bipartite Gaussian state. We derive a pre- 
Lindblad master equation in the non-rotating wave approximation for the system. We then solve 
the master equation for a bipartite harmonic oscillator Hamiltonian with entangled initial state. We 
show that for strong damping the loss of entanglement is the same as for freely evolving particles. 
However, if the damping is small, the entanglement is shown to oscillate and eventually tend to a 
constant nonzero value. 
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Motivation 

Entanglement is one of quantum mechanics' most 
fascinating features. It was first described in a celebrated 
paper by Einstein, Podolsky and Rosen [l[ but owes its 
name to Schrodinger [2[ , who investigated its broader sig- 
nificance for the measurement question. It has taken on 
enhanced significance in quantum information. In this 
regard, the fragility of entanglement when the system is 
subjected to "outside" influence is of even greater impor- 
tance. In the current work, we study a bipartite system 
with a Gaussian wave function. The state is prepared 
such that it is entangled, then shared between two par- 
ties who let their respective particle evolve either freely 
or interacting via a harmonic potential, but interacting 
with its own environment or heat bath. We study the 
resulting loss of entanglement between the particles. To 
do so, we first give a simple derivation of a pre-Lindblad 
non-rotating- wave master equation, 0, Q , starting with 
the Quantum Langevin Equation as derived in |5| and 
using a simple perturbation method as in Q. 

The loss of entanglement in a system interacting with 
an environment is a well-known phenomenon. It has 
been studied in various systems, see e.g. 043' where 
it was found that there is often a sharp loss of entangle- 
ment when compared to a decoherence time scale, which 
has been termed entanglement sudden-death (E.S.D.). 
These studies are mainly in the context of qubits and 
the Rotating Wave Approximation (R.W.A.) however, 
whereas this work presents a study of E.S.D. in a 
continuous-variables setting and uses the Non-Rotating- 
Wave (N.R.W) approximation. Note that the master 
equation obtained in the N.R.W approximation is not 
of the Lindblad form [l3| , hence could in principle have 
unphysical results. However, the unphysical behavior is 
obtained for low T, and one can easily check for the va- 
lidity of the density matrix by checking its positive semi- 
definiteness. On the other hand, the N.R.W master equa- 
tion often works better for systems which are strongly 
coupled to the environment [3(. 



In [fjj , Ficek and Tanas study a system of two qubits 
coupled to a radiation field where they allow sponta- 
neous decay of the atoms. They show that the entan- 
glement vanishes but then is revived twice. In [l5j . the 
authors study the emergence of entanglement between 
two initially non-entangled qubits due to spontaneous 
emission, provided both atoms are initially excited and 
in the asymmetric state. Their results suggest that an 
interaction between two particles which are initially en- 
tangled can delay the vanishing of the entanglement and 
even revive it, or create entanglement between two ini- 
tially non-entangled particles. We introduce a harmonic 
potential with frequency ujq as the interaction between 
the particles in our system and examine the dynamics of 
the entanglement. We show that entanglement revival 
can occur depending on the strength of the damping, i.e. 
how strong the coupling 7 is with respect to the oscil- 
lator's frequency. We show that if the damping is small 
(7 < 2\/2ljo), the entanglement oscillates towards a cen- 
tral value and does not vanish asymptotically. 

In Section U we recall the Langevin equation and the 
main steps in the derivation of the master equation. We 
then establish in Section [XT] the formalism used to de- 
scribe Gaussian states and the particular measure for 
entanglement we use. Section 11111 illustrates E.S.D. while 
Section HVl contains the main results of this paper. Sec- 
tion [V] contains some concluding remarks. 



FRAMEWORK 



The derivation of the master equation is detailed in 
more detail in Appendix [V] and will only be succintly 
presented in the following. The derivation is easiest for 
one-particle but generalises just as easily to the case of 
two-particles, each coupled to its own environment. We 
consider a heat bath modelled by independent oscillators 
coupled harmonically to the particle [5| . The correspond- 
ing Hamiltonian has the form 
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Solving the Heisenberg equations of motion for qj yields 
the Quantum Langevin Equation 

mx + I p(t - t')x(t') dt' + V'(x) = £(t) (2) 

J — oo 

where the dot denotes the derivative with respect to time 
and the prime that with respect to x. p(t) and £(t) de- 
scribe the influence of the bath on the system and are 
known as the memory function and the operator- valued 
random force respectively and are expressed explicitly in 
Appendix |Vj In the case of a Ohmic heat bath, p,(t) 
effectively reduces to 7 and we can write, for a general 
observable Y of the small system (particle) , the Quantum 
Langevin Equation reads 

Y = ~ [H s , Y] - ± [[as, Y] , £(t)] + + ^ [[x, Y] , x(t)]+ ■ 

(3) 

For a system of two particles, each connected to their 
individual heat bath, we have analogously 

Y=\[H S ,Y ] 

-^[[^Y]Mt)] + -^[[^Y]Mt)]+ 
+ 2k Y] ' 7lil(<)] + + ih [[a;2 ' r] ' ^2(t)}+ 

(4) 

The Langevin equation is an equation for the system 
operators (Heisenberg representation), whereas a mas- 
ter equation is an approximate equation acting on the 
density operator of the quantum system under study 
(Schrodinger picture). The adjoint equation provides a 
link between the two formalisms. We define 

Tr s {Y(t)p} = Tr s {Y pit)} 

and obtain the adjoint equation 

Ps{t) = - l -[H s , Ps (t)]-^_[[i 1 {t), Ps {t)] + ,x 1 ] 

~[Mt),Ps(t)] + ,x 2 ] (5) 

+ 2^ [[±l,Ps(t)] + ,Xl] + ^ [[±2,Ps(t)} + ,X 2 ] 

In order to derive the master equation let us assume that 
the bath is large so that we may assume that it stays 
at thermal equilibrium and that at t — > —00, the sys- 
tem and the bath are decoupled so that po(t) = p s {t)pB- 
This assumption is critical to the derivation of any mas- 
ter equation. Finally, assuming that the noise is small 
allows us to write — > e£(t). This assumption is not 
essential to the derivation but allows for a simpler deriva- 
tion. Applying a perturbation method and tracing over 
the bath yields the Non-Rotating- Wave master equation 



for p s (t) 

m = -\{Hs,p{t)\ 

+ 7^ [[xi,p(t)] + , xt] + -^r [[± 2 , p(t)} + , x 2 ] 

(6) 

II. GAUSSIAN STATES AND THE 
LOGARITHMIC NEGATIVITY 

Since the states we will study are Gaussian, we now 
briefly recall the formalism for Gaussian states [l6l - [l8l |. 

Gaussian states can be completely specified in terms 
of their first and second moments, described respectively 
by the displacement vector 

d 3 = ( J R i ) p = Tr[R jP ] 

and the covariance matrix 

7j , k = 29^^ - (Rj) p )(R k - (R k ) p )} 

where R is the vector R T = (gi, pi ;....; q n ,Pn) ; <lj and 
Pj are the canonical variables of a system of n oscillators 
with the usual canonical relations written as [Rj,R k ] = 
ihajk and a a real skew-symmetric 2n x In block matrix 
given by 

k=l v 7 

The displacement vector are irrelevant in the study of 
entanglement and are taken to be zero in our examples. 
The covariance matrix thus reduces to 

■2Kc\vpR J lh (7) 

Any real symmetric positive-definite matrix A can be 
brought to its Williamson normal form [f9j via symplec- 
tic transformations, i.e. transformations that preserve 
the canonical commutation relations, Awf = SAS T = 
diag(ai, a±, ....a„, a n ) where the a^'s are the symplectic 
eigenvalues of A. One can calculate them as the positive 
eigenvalues of iaA or more simply as the positive square 
root of the eigenvalues of —crAcrA. 

A particularly suitable measure of the entanglement of 
mixed Gaussian states is the logarithmic negativity 
HEEIS- It vanishes for separable states, does not increase 
under LOCC (local operations and classical communica- 
tion), and stays invariant under local unitary transfor- 
mations. It is defined as 

In 

£^(p) = -]Tlog 2 (mm(I,| A, I)), (8) 
1=1 
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where the \ are the symplectic eigenvalues of the par- 
tially transposed covariance ^ Tl \ which is obtained from 
7 by reversing the time in all variables of one of the sub- 
systems. Choosing to transpose with respect to particle 
1, we replace x\ — > x\ and p\ — > —p\. The A^'s are thus 
the square roots of the eigenvalues of —aj Tl aj Tl . 



which, when substituted into (fTOI) yields 



P(q,z,t) 



-A 1 q 1 2 -A 2 q 2 -BiZ-i z -B 2 z 2 z -Dz 1 z 2 -Eq 1 q 2 



-Cuz 1 q 1 —C 22 z 2 q 2 —Ci 2 z 1 q 2 -C 2 iz 2 qi 



(14) 



where the coefficients are given by 



III. FREE EVOLUTION OF AN ENTANGLED 
INITIAL STATE 



We first consider the case of a free particle Hamitonian 

H = ?L + A (9) 

2m 2m v ; 

with p\x) = —\h-g^\x). This will allow us to examine the 
dynamics of the entanglement when an entangled bipar- 
tite Gaussian state is left to evolve, each particle coupled 
to its own heat bath. ([6]) can be solved to obtain 

P(q, z, t) = P(q, z (i), 0) exp [-75 q\ t - n q\ t] 
xexp -Ai(t) + -^(*) (** + ^ 

+01 (t) [ zi + ^- ) + a 2 (i) ( zi + 



with 
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Ai(t) = 2mfcT 4 (l - e - 2 ^*/™) 
a .( t ) = 4sm.(i_ e -w*/tn) 



272 



(10) 



7» 

r,: 



zo,i(*) = z» e 



kTj 

7i 

27; I 



-7*t/r. 



The full derivation for this solution can be found in 
Appendix |Vj Let us consider a bipartite initial state 
with the Gaussian wavefunction, suggested by Ford and 
O'Connell ED, El 



\&( Xl ,x 2 ) = f2 1/2 e" <X1 4^ 2> e 
The corresponding density matrix is 



(»i+«»r 



(11) 



p| = f^g- e +( :r i 2 + :E 2 2 +a;' 1 2 +X2 2 )+2e_ (s 1 a2+a:iX2) ^2) 



where ft = ^ e± = ^ ± 



We next compute the time-evolved state by insert- 
ing the Fourier transform of (TT21 into (flU)) . The Fourier 
transform of (fl"2l is 



P(q z ;0) 

= exp [— 2e + h 2 zoi — 2e + h 2 zo\ + 4e_ft 2 Z01Z02] 



x exp 



e+(<Zi 2 + ?2 2 ) 



e-<7i<?2 



(e^-el) 4(e 2 h -e 2 _) 



(13) 



eP s 2 ctj Xj h 2 e+ 



27 . ■ 4 7 .2 ■ 27 . 2 (1 e 



P, = 2fi 2 e+e" 1 ^- + Xj 

_ 2S, 2 e ± _7£6 _2t£s /A., . 

7j V7j ' 



with 



e+ ifj = k 
e_ if j 7^ fc 



P = 2?i 2 e_e 



P = 2d 2 



Jit _ 12 t 



f -71 t 79 1 

£ --(l-e-^)(l-e-^) 



7i 72 



(15) 



The entries of the covariance matrix can be calcu- 
lated directly from (fl"4]) taking into account the change 
of variables (l44l): 



-KI^ (, ' z=0 ' t,l| •■- 

d d ~ 
— — P(q,z,t)| q=0 , z= 

1/5 a 



2Re(X i X j ) 
2Re(X i P i ) 
2Re(P i A j ) 
2Re(pP i ) 



2 V dzj dzj 



P(q=0,z,i) | z=0 



The covariance matrix is then 



It = 



( 


4Ai 


-C11 


E 


— C21 


-C11 


Pi 


— C\2 


D 




P 




4A 2 


—C22 


\ 


— C21 


D 


— C22 


P 2 



(16) 



We now perform the partial transposition with respect 
to particle 1: 



which is real and symmetric. 

The symplectic eigenvalues are given by 

A ± = ± - V( £ ll - £ 33j 2 + 4£i 3 £24 - 4ei 4 e 2 3 



4Ai 


C n 


E 


— C21 


C11 


Pi 


C\2 


-P 


P 


C12 


4A 2 


— C22 


— C'21 


-P 


—C22 


P 2 



(18) 
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Logarithmic Negativity for free evolution vs time 



gm1 = 


1,gm2 = 


1 , s= 0.25 


gm1 = 


1, gm2 = 


1,8=1 


gm1 = 


1, gm2 = 


1,8=2 



FIG. 1: Logarithmic negativity vs t for three values of s 

The values of s are : dashed s = 0.25, dotted s — 1, 
dash-dotted s = 2 



where E\ 

Ell = 
£33 = 
£13 = 
£14 = 

£23 = 
£24 = 



£21 = £34 = £43 = and 



£22 
£44 
£42 
-£32 
'£41 
£31 



AAyBi -DE + C 12 C 21 - C\ x 
AA 2 B 2 ~Cl 2 - DE + C 12 C 21 
EB\ — 4A 2 D — C\\C\ 2 + C\ 2 C 22 
—C\ 2 B 2 — C 2 \B\ + C\\D + C 22 D 
-EC n + iAtOa + ^A 2 C 21 - EC 22 
EB 2 — C 22 C 2 \ + C\\C 2 \ — 4AiD 



(19) 



The logarithmic negativity then becomes 

£}<{p) = -2 (log a (min(l, |A^|)) + log 2 (min(l, |A^|))) 

(20) 

Figure Q] shows the logarithmic negativity as a function 
of time for three values of s. We can observe that there 
is complete disentanglement between the particles from 
a sharp cut-off time onwards, which obviously depends 
on s, and hence on the initial degree of entanglement. 
The sharp cut-off time characterizes entanglement sud- 
den death (ESD). 



IV. EVOLUTION WITH A HARMONIC 
POTENTIAL INTERACTION 



pendixfV] In this case 
/ 



M = 
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— 4m 2 WQ 
V 4m 2 w 2 




4m 2 WQ 

-4m 2 w 2 



(22) 



In general, the eigenvalue equation is quartic and the 
solution is complicated. We therefore assume, for sim- 
plicity, that 71 = 72 = 7 and T\ = T 2 = T. In that case 
the eigenvalues are 



^0, 27, 7 + yff - 8™ 2 w 2 n-\ff 



(23) 



and after some unpleasant algebra, we can write 

P = exjp[-Aql - Aql - Eqiq 2 - Bz\ - Bz\ - Dz\z 2 
-Cizxqx - Ciz 2 q 2 - C 2 zxq 2 - C 2 z 2 qi] 

(24) 

which is of the same form as (fl"4"]l . with A\ = A 2 = A, 
B 1= B 2 = B, C11 = C 22 = Ci and C 12 = C 21 = C 2 ex- 
cept that the explicit expressions for A, B, etc. are more 
complicated, and will be omitted here. The following fig- 
ures show the logarithmic negativity for various values of 
wo and of 7. 

Logarithmic negativity with and without ootential in the over-damped case 



3 = 3,0=1 
gm1=3 : gm2 = 3 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 2: Logarithmic Negativities with and without potential 
for uq — 1 and 7 = 71 = 72 = 3 

The dotted line represents £ plotted with the potential. 



If we introduce a harmonic potential interaction into 
the Hamiltonian, ([9]) generalises to 

„ Pi , Pi , ™Ai N2 foU 

Hs = 2^ + 2^ + — {Xl - X2) (21) 

We can include this into © and solve the resulting dif- 
ferential equation following the method described in ap- 



V. CONCLUDING OBSERVATIONS 

Figurc[2]and Figure[3]illustrate that, in the presence 
of a harmonic interaction between the particles, there is 
a marked difference in behaviour between two damping 
regimes. In the over-damped case (7 > 2V2wo), Fig- 
ure [U it can easily be seen that the curves coincide. This 
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Logarithmic negativity in the under-damped case 



g= 1,o= 1.5 

g= 1,o= 1.8 
g= 1,0=2.2 




FIG. 3: Logarithmic Negativities with and without potential 
for ujq = 1 and 7 = 71 = 72 = 0.2 

The dotted line represents 8 plotted with the potential. 



Logarithmic Negativity with the harmonic potential 



Log Neg 


3=1,0=1.5 


Log Neg 


3= 1,o= 1.8 


Log Neg 


1,o- 2.2 




FIG. 4: £ in the highly under-damped case 

The plots are obtained with 7=1 and : full ujq — 1.5, 
dashed ujq — 1.8 and dotted ujq = 2.2 



suggests that if the coupling is much stronger than the 
harmonic potential, the decay of the entanglement is un- 
affected by the potential. As a matter of fact, closer 
investigation of the vanishing point shows that the en- 
tanglement disappears faster as damping increasing. On 
the other hand, if the damping is small (7 < 2^/2ujo), the 
entanglement can reappear several times, as illustrated 
by Figure [3] Figure [4] shows the behaviour of the entan- 
glement as the damping decreases still. It can easily be 
seen that as the harmonic potential becomes stronger, the 
entanglement does not disappear. Instead it decreases 
sharply before being "restored". It then tends towards 
a non-zero constant value for large times. This suggests 
that allowing the particles to interact harmonically effec- 
tively saves the entanglement. 
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The Master equation 

In this section, we will describe the derivation of the 
master equation in more details. A simple derivation of 
the Quantum Langevin Equation, starting from a heat 
bath modelled by independent oscillators coupled har- 
monically to a system of one particle, was given in Q. 
The Hamiltonian has the form 



H 



P 

2m 



^-+ m ^(q J -x) 2 \ (25) 



Solving the Heisenberg equations of motion yields 



Qi{t) = cos(w,i) -I sin(cj ? i). 

Introducing the quantities 

li(t) = ^ m jU 2 j cos(ojjt)Q(t) 

3 



cos [u>j(t -t')]x(t')d£ 

(26) 



(27) 



(0(i) is the Heaviside function) we obtain the Quantum 
Langevin Equation 

mx k + I n k (t - t')x k (t') dt' + V'(x k ) = £fc(*X28) 

J —CO 

where the dot denotes the derivative with respect to time 
and the prime that with respect to x. fi(t) and £(t) de- 
scribe the influence of the bath on the system and are 
known as the memory function and the operator-valued 
random force respectively. We also introduce the spectral 
distribution 

G{uj) = mc [fi(uj + iO+)] 

- ^Y^m^lSiu-u^ + Siu + Ui)], (29) 
3 

in terms of which the autocorrelation of £ (t) is given by 
1 f°° f Too \ 

- / G{uj)huj coth I — — ) cos [uj{t - t')] duj (30) 
7T Jo \2kTJ 
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where [, ]+ denotes the anticommutator. For a general 
observable Y of the small system (particle) , one can write 



y = - n [H,Y] 



UH.,Y]-±[[x,Y],m + 



1 



2h 

[x,Y] , / d£ fi(t - t')x(t') 



(31) 



In the case of an Ohmic heat bath, we can replace 



p(t')x(t') dt' -> 7±(t) and G(u) -> 7 (32) 



so that the Quantum Langevin Equation reads 



Y = - h [H s ,Y] 



-[[x,Y],mU + ^[[^Y],x(t)] + 



2h 



2h 



(33) 

For a system of two particles, each connected to their 
individual heat bath, we have analogously 



Y 



:[H S ,Y] 

— y] , ex (t)]+ - 4 t^, r i ^ 2 (*)]+ 



2ft 



+ ^ [[H, >1 , 71*1 («)]+ + ^ [[S2, ^ , 72*2(«)] 



+ 

(34) 



If we define 



can be solved successively. The equation for po reads 



p (i) = --[H.,po(t)] 



+ 2^ [[xi,Po(t)] + ,xi] [[x 2 ,po(*)] + ,a!2 



2fi 



(37) 



The equation for p! can be written as 



p x (t) = --[H s , Pl (t)] 



+ 2E H^P^h^ 1 ] - 2ftKi(*)'Ps]+[po(t),a;i] 

+ 2§ [[*2, pi (*)]+, £2] - 2^[6(*),PB]+[po(t),a;2]. 

(38) 



The solution for the first-order term can be written as 
rt 

>PBj + 



PiW = ~J e A ^{ {p (t'), Xl ] ® 

+ |>o(*').a!2]®[6(0»PB] + }<tt / 

where v4 s is a super operator which, applied to pfc(i) 
yields 



A s p fe (t) = --[ff s ,p fe (t)] 



+ 2^ [[*i,Pfc(t)] + ,a;i] + ^| [[i 2 ,Pfc(i)]+ ,^2] 



2/) 



(39) 



Finally, we insert this solution into the equation for p 2 : 



Tr s {Y(t)p} = Tr s {Yp(t)} 



(35) 



where Tr s is the trace over the system. Let us introduce 
p(t) = p s {t) ® ps where p s is the density matrix of the 
system and ps that of the bath. It follows easily from 
that p(t) satisfies the adjoint equation 



^[^..pW]-^[Ki(*).P(*)]+.a:i] 
k [Mt),p(t)}+,X2] + 1 ^-[[i 1 ,p(t)} + ,x 1 ] 



P(t) 



2h 

— [[x 2 ,p(t)} + ,x 2 ] 



(36) 



To derive the master equation (that is, an effective equa- 
tion for p s {t)) from the adjoint equation, we assume that 
the noise is small and temporarily introduce a small pa- 
rameter e, replacing by e£(t). We can write v{t) to 
second order in e as 



hit) = --[H s , P2 (t)\ 



+ zt [[^'^WU' 2 ^] +"dr [[*2,P2(t)] + ,a;2 



2h 



- 2ft [[&(*). Pi (*)]+> *i] - ^ [fe (*) ' Pi (*)]+ > *a] 

(40) 

Taking the trace over the bath variables and using the 
autocorrelation (|30[) in the Ohmic limit, we obtain the 
non-rotating-wave master equation (upon removal of e) 



Ps (t) = --[H s ,p s (t)] 



171 
' 2h 
fcTVyi 
h 2 



+ 7^ [[xi,Ps(t)] + ,xi] + ^| [[ar 2 ,p a (t)] + ,x 2 ] 
[p a (t),xi] ,aci] - fc ^ 2 2 72 [[p s (t),x 2 ] ,x 2 ] 



fi 2 



(41) 



We also assume the baths and the system are decoupled 
at t = —00, so that po(t) — po(i)ps- Inserting this ex- 
pansion into (|3l)]) yields equations for po, pi and p 2 which 



Solution to the master equation 

The derivation of the solution to the master equation 
will here be given in a general way. The method will be 
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described for a free particle system Hamiltonian 

H= A + A (42 ) 
2m 2m v ; 

with p\x) = —\h-^\x), but is easily generalised to other 
types of Hamiltonians. Note that we assume that the 
particles have the same mass. In position-space, © be- 
comes 



d_ 

at 



'Mp\v) = 

m ( d 2 d 2 
2m \dx 2 dy\ 



dx\ 



«~2 I P 



{xi -yi) p p— {x2-y%) p 

Using the change of variables 

x = u + hz and x = u — hz 



(43) 



(44) 



and replacing p(x, x', 0) — > P(u, z, 0), we apply a Fourier 
transformation with respect to u: 

P(q,z,f) = Ar [ P(u,z,t)e~ iqiUl - iq2U2 du 1 du 2 (45) 
4tH J 

obtaining an equation for P(q, z, t): 



v rn 

"72 



d_ 

Q2_\ _d_ 



- Zl+ 2^J^ +47lfcTlZl 



: ^2^J^ +472fc ^ 



P(q,M) (46) 
P(q,z,t).(47) 



This equation can in principle again be solved using the 
method of characteristics. The characteristic equation is 



dv _ M 
~dt ~ 2m 

with v = Z2, qi, q%) T and 



(48) 



M 



On a characteristic, 
^P(q, *(*),*) 



/ 2 7 i 10 

2 7 2 1 



V 



(49) 



= -[^mzlit) + A l2 kT 2 z 2 2 {t)]P{ci,z{t),t) 



The eigenvalues and eigenvectors of M can be computed 
to be 



and 



A T = (2 7 i,27 2 , 0, 0) = (Ax, A 2 , A 3 , A 4 ) (50) 



/I 


1 

271 


^ 


1 





1 

272 





1 





\ 





1 / 



(51) 



Since Q 1 MQ = D where D is the diagonal matrix, we 
need Q~ x as 



Q- 



( l wA 

1 4 

1 

V )o o o i / 



(52) 



Then we can write 2m^ = flw with w = Q x v which 
is easily solved so that v(t) = Qe Dt ' 2m Q~ 1 Vo with 



,Dt/2m 






V o 



o o\ 

e J2t/m o o 



/ 



(53) 



Some more algebra yields the solution 

P(q, z, t) = P(q, z (t), 0) exp [-r 2 g 2 t - n g 2 t] 
xexp -Ai(t) + ( Z2 + || 

+ax(t) + +«>(*) (,x + ^ 



(54) 



with 



Ai(t) = 2mkT.il - e - 2 7i*/m) 
a .(i) = 4mMl(i_ e -7.t/™) 



«b,<(t) = «i e- 7 ' t/m - (1 - e-^/ m ) 



We thus have an expression giving the time dependency 
for an arbitrary initial state. 
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